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O ■ ABSTRACT 

^ . HH 1 10 is a rather peculiar Herbig-Haro object in Orion that originates due to the deflec- 

^ ■ tion of another jet (HH 270) by a dense molecular clump, instead of being directly ejected 

from a young stellar object. Here we present new results on the kinematics and physical con- 
ditions of HH 110 based on Integral Field Spectroscopy. The 3D spectral data cover the whole 
^ ■ outflow extent (~ 4.5 arcmin, 0.6 pc at a distance of 460 pc) in the spectral range 6500-7000 

Cn\ A. We built emission-line intensity maps of Ha, [N ll] and [S ll] and of their radial velocity 

(N . channels. Furthermore, we analysed the spatial distribution of the excitation and electron den- 

sity from [Nll]/Ha, [S ll]/Ha, and [S ll]67 16/6731 integrated line-ratio maps, as well as their 
behaviour as a function of velocity, from line-ratio channel maps. Our results fully reproduce 
■ the morphology and kinematics obtained from previous imaging and long-slit data. In addi- 

I tion, the IFS data revealed, for the first time, the complex spatial distribution of the physical 

. conditions (excitation and density) in the whole jet, and their behaviour as a function of the 

kinematics. The results here derived give further support to the more recent model simulations 
that involve deflection of a pulsed jet propagating in an inhomogeneous ambient medium. The 
. , IFS data give richer information than that provided by current model simulations or laboratory 

/\ I jet experiments. Hence, they could provide valuable clues to constrain the space parameters 

^ . in future theoretical works. 
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1 INTRODUCTION 

HH 110 is a Herbig-Haro (HH) jet emerging from the southern 
edge of the L1617 dark cloud in the Orion B complex. Both ob- 
servational a nd theoretical works have been carried out since its 
discovery by iReipurth & OlbergI (Il99lh , because unlike most of 
the well-known stellar jets, HH 110 has a peculiar morphology, 
among other properties. The morphology of HH 1 10 in the Ha and 
[S II] lines is very complex, starting in a collimated chain of knots. 



E-mail: rosario@am.ub.es; bgarcia@iac.es ; sanchez@cefca.es; 
gabriel . gomez @ gtc . iac . es ; robert . estalella @ am . ub. es ; angels .riera @ upc . edu 
t Based on observations collected at the Centro Astronomico Hispano 
Aleman (CAHA) at Calar Alto, operated jointly by the Max-Planck Institut 
fiir Astronomie and the Instituto de Astrofisica de Andalucia (CSIC) 



The emission away from the collimated region has a more chaotic 
structure and widens within a cone of unusually large opening an- 
gle 10°). HH 110 appreciable wiggles along the ^ 4 arcmin jet 
length. The knots are all embedded in fainter emitting gas, which 
outlines the whole flow, more reminiscent of a turbulent outflow. 
Most of the knots detected in ground-based images are spatially 
resolved into several components in the h igher spatial resolution 
HST images (see e. g. iHartigan et al.ll200^ . The transverse cross- 
section of HH 110 shows a significant asymmetry, the eastern bor- 
der is sharp and poorly resolved, whereas the strong knotty emis- 
sion mostly appears towards the western side. HH 110 is the only 
known HH that shows a faint filament of emission lying parallel at 
^10 arcsec to the east of the northern A-C knots, which prob- 
ably represents a weak secondary or even a fossil flow channel 
CReipurth & Olberg. . 1 99 1.) . Attempts to find the driving source of 
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HH 110 have failed at optical, near infrared and radio continuum 
wavelengths. 

A scenario that accou nts for the singular mor phology of 
HH 110 was first outlined bv lReipurth. Rasa & Heathc ote (199^). 
These authors proposed that HH 110 originates from the deflec- 
tion (deflection angle ^ 60°) of the adjacent HH 270 jet through a 
grazing collision with a dense molecular clump of gas. The feasi- 
bility of this scenario has been further reinforced from the results 
of numerical simulations, which model the emission a rising from 
the co llision of a jet with a dense molecular clump (see lRaga et al.l 
I2OO2I and references therein). In addition, analysis of further data 
are also in agreement with this scenario. First, a high-density clump 
of gas around the region where the HH 270 jet changes its di- 
rection to emerge as HH 110 has bee n detected through high- 
de nsity tracer molecules (in HCO+, bv IChoil[200ll . and in NH3, 
bvlSepulveda et al.ll201Qh . Second, results from proper motion de- 
termination are also consistent with t he jet/clo ud collision scenario 
jReTpurth. Raga & Heathcotd 1 19961 : iLopez et al.ll2005h . A recent 
work of lHartigan et al. (200^) that included laboratory experiments 
has given further support to the jet/cloud collision scenario. 

HH 110 is a good candidate to search for the observational 
footprints of gas entrainment and turbulence by analysing the kine- 
matics and the excitation conditions along and across the jet flow. 
Some works were performed in the recent past from long-slit spec- 
troscopy and Fabry-Perot data. The kinematics and physical condi- 
tions, both along the outflow axis and at four positions across the 

i' et beam , were explored from long-slit spectro scopy by Riera et all 
2003ah . who found very complex structures. iRiera et all (l2003bh 
explored the spatial distribution and the characteristic knot sizes, 
as well as the spatial behaviour of the velocity and line width, by 
performing a wavelet analysis of Fabry-Perot data, but only cover- 
ing the Ha line. Their results indicated that most of the Ha kine- 
matics can be explained by assuming an axially peaked mean flow 
velocity, on which are superposed low-amplitude turbulent veloci- 
ties. In addition, their results are suggestive of the presence on an 
outer envelope that appear s to be a turbulent boundary layer. Fi- 
nally, iH^tiganetal] J2009h compared images from laboratory jet 
experiments with numerical simulations and with long- slit, high- 
resolution optical spectra obtained along HH 110. They found a 
good agreement between the shock structures observed in HH 110 
and those derived from experiments of a supersonic jet deflected by 
a dense obstacle. 

In order to further advance on our understanding of HH 110, 
this object was included as a target within a program of Inte- 
gral Field Spectroscopy of Herbig-Haro objects using the Potsdam 
Multi- Aperture Spectrophotometer (PMAS) in the wide-field IFU 
mode PMAS fibre PAcK (PPAK). The data obtained from this IFS 
HH observing program give a full spatial coverage of the HH 110 
emission in several lines (Ha, [Nil] and [Sll]), thus allowing us 
to perform a more complete analysis of the kinematics and physi- 
cal conditions through the whole flow than all the previous works. 
The main results are given in this work. The paper is organized as 
follows. The observations and data reduction are described in § 2. 
Results are given in § 3: the analysis of the physical conditions in § 
3.1 and the analysis of the kinematics in §3.2. A summary with the 
main conclusions is given § 4. 



2 OBSERVATIONS AND DATA REDUCTION 

Observations of HH 110 were made on 22 November 2004 with 
the 3.5-m telescope of the Calar Alto Observatory (CAHA). Data 




Figure 1. CCD image of HH 110 obtained at the Nordic Optical Telescope 
(NOT) through a narrow-band filter that includes the [S ii] 6716, 6731 A 
lines (see Lopez et al. 2005 for details). The four IFU pointing fields are 
superposed on the image. The knots discussed throughout this work are 
labeled on the image, ac cording to the nomenclature first established by 
iReipurth & 01berdll99lh . The contours of the individual fibres have been 
drawn on the southern IFU field. 



were acquired with the Integral Field Instrument Potsd am Multi- 
Aperture Spectrophotometer PMAS jRoth et all l2005l) using the 
PPAK configuration that has 331 science fibres, covering an hexag- 
onal FOV of 74x65 arcsec^ with a spatial sampling of 2.7 arcsec 
per fibre, and 36 additional fibres to sample the sky (see Fig. 5 in 
iKelz et al.ll200^ . The 11200 grating was used, giving an effective 
sampling of 0.3 A pix~^ (^ 15 km s~^ for Ha) and covering the 
wavelength range ^ 6500-7000 A, thus including characteristic 
HH emission lines in this wavelength range (Ha, [N II] AA6548, 
6584 A and [Sll] AA6716, 6731 A). The spectral resolution (i.e. 
instrumental profile) is ^ 2 A FWHM (^ 90 km s~^) and the ac- 
curacy in the determination of the position of the line centroid is 
'-^ 0.2 A (^ 10 km s~^ for the strong observed emission lines). 



Table 1. Positions of the single pointings of the HH 110 mosaic. 
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Position 


Exp. time(s) 


knot 


05^51^25^3 


+02°55'16'' 


1800 


A-E 


05h51™24!4 


+02°54'20'' 


1800 


I-N 


05h51™23!8 


+02°53'41'' 


2400 


N-R 


05h51"^24!9 


+02°52'40'' 


2700 


S 



Four overlapped pointings were observed to obtain a mosaic of ^ 
5^x \'5 to cover the entire emission from HH 110 (see Fig. [T]). 
Table 1 lists the centre positions of each pointing, the exposure 
time and the HH 1 10 knots included in each po inting, following 
the nomenclature o f iReipurth. Raga & Heathcot3 Jl996h . 

Data reduction was performed using a preliminary version of 
the R3D software ( Sanchez 2006), in combination with IRAlQ and 
the Euro3D packages (Sanchez 2004). The reduction consists of the 
standard steps for fibre-based integral field spectroscopy. A mas- 
ter bias frame was created by averaging all the bias frames ob- 
served during the night and subtracted from the science frames. 
The location of the spectra on the CCD was determined using 
a continuum-illuminated exposure taken before the science expo- 
sures. Each spectrum was extracted from the science frames by 
co-adding the flux within an aperture of 5 pixels along the cross- 
dispersion axis for each pixel in the dispersion axis, and stored in a 
row- stacked- spectrum (RSS) file (Sanchez 2004). The wavelength 
calibration was performed by using the sky emission lines found 
in the observed wavelength range. The accuracy achieved for the 
wavelength calibration was better than ^ 0.1 A 5 km s~^). 
Furthermore, the large scale diffuse emission from this Orion re- 
gion was subtracted from our data by using the signal acquired 
through the 36 additional sky fibres. Observations of a standard 
star were used to perform a relative flux calibration. The four IFU 
pointings were merged into a mosaic using our ow n routines, de- 
veloped for this task (see e. g. Sanchez et al. 2007 and references 
therein). The procedure is based on the comparison and scaling 
of the relative intensity in a certain wavelength range for spatially 
coincident spectra. The pointing precision is better than 0.2 arc- 
sec, according to the pointing accuracy of the telescope using the 
Guiding System of PMAS in relative offset mode (used for these 
observations). The overlap between adjacent pointings is ^ 10- 
15% the field of view (more than 30 individual spectra). The merg- 
ing process was checked to have little effect on the accuracy of 
the wavelength calibration of the final datacube making use of the 
sky lines (i. e. a common reference system) present in the spec- 
tra. A final datacube containing the 2D spatial plus the spectral in- 
formation of HH 110 was then created from the 3D data by us- 
ing Euro3D tasks to interpolate the data spatially until reaching 
a final grid of 2 arcsec of spatial sampling and with a spectral 
sampling of 0.3 A . Further manipulation of this datacube, de- 
voted to obtain integrated emission line maps, channel maps and 
position-velocity maps, were made using several common-users 
tasks of STARLINK, IRAF and GILDAS astronomical packages and 
IDA (^rcia-Lorenzo, Acosta-Pulido & Megias- Fernandez. .2002.) . 
a specific IDL software to analyse 3D data. 



IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 



3.1 Physical conditions 

3.1.1 Morphology 

Some morphological differences have been found in previous 
narrow-band images of HH 110 between the Ha-i-[Nll] and [S ll] 
emissions. In an attempt to explore these differences, we obtained 
IFS-derived narrow-band images for each of the emission lines in- 
cluded in the observed spectral range (Ha, [Nil] A 6584 A and 
[Sll] A 6716,6731 A). For each position, the flux of the line was 
obtained by integrating the signal over the wavelength range of the 
line (A 6558.3-6565.9 A for Ua; A 6579.6-6585.8 A for [Nil]; A 
6712.8-6718.1 A for [Sll] A67 16 and A 6727.2-6733.1 A for [S ll] 
A6731), and subtracting a continuum, obtained from the adjacent 
wavelength range free of line emission. The narrow-band maps of 
HH 110 obtained from the IFS data (Fig.|2]) are in good agreement 
with those obtained for this jet through gr ound-based, narrow-band 
imag es by different authors (see e. g. Reipurth, Rag a & Heathcotd 
1996; Lopez et al. 2005). However, since the IFS data allow us to 
properly isolate the Ha line emission from the [N II] lines, we ob- 
tained the first known image of HH 110 in the [Nil] emission. In 
this sense, the IFS maps are better suited for comparing the jet mor- 
phology in the Ha and [S II] emissions, since the Ha emission is 
not affected by the contamination from [N II] lines. 

Most of the HH 110 knots were detected in all the emission 
lines, although the brightness varies among different lines, being 
a signature of the excitation conditions, as we will discuss later. 
There are a few emission features whose brightness in the [Nil] 
and [S II] lines should be very weak relative to the Ha brightness 
and were barely detected or not detected at all in the [Nil] and 
[Sll] lines, namely: (i) the knot K, appearing quite isolated to- 
wards the western outflow edge, but outside the lower-brightness 
emission surrounding the knots, and (ii) the extended features to- 
wards the east of knots A to C, labeled X and Y. These faint Ha 
featur es were first detected in th e ground-based, narrow-band im- 
age of Reip urth & Olberd (Il99lh . Our new IFS data indicate the X 
and Y features to be a "true" emission in the Ha line, having nei- 
ther significant contributions from the extended nearby continuum 
(that was removed) nor from [N ll] lines (if any, it should be below 
the 3(7 rms level). 

The knots are found surrounded by a lower-brightness nebular 
emission for all the emission lines. This more diffuse emission also 
shows a shock-excited spectrum, and likely arises from the gas in 
the boundary layer that is being incorporated into the outflow. The 
large scale, low-velocity diffuse emission component detected in 
the Orion region was removed from our data and is not contribut- 
ing to this low-brightness, interknot emission structure we are re- 
ferring to. Figured shows that the cross section along all of the jet 
appears wider in Ha than in the other lines, which all have similar 
spatial widths. Since with the present spatial resolution the sizes of 
the knots are similar at all lines, the wider Ha cross-section of the 
jet should arise from the contribution of the low-brightness, nebu- 
lar emission. In our maps, this weaker emission appears spatially 
more extended in Ha than in the rest of the observed lines. This 
fact could be due, in part, to the weakness of the [N ll] and [S ll] 
emissions, which makes the emission undetectable in our maps (up 
to the 3(7 rms level) farther away from the knots. 
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Figure 2. Integrated maps of the HH 110 emission obtained from the datacube by integrating the signal within the wavelength range including the line labeled 
in each panel: (Ha : A6558.3-6565.9 A, [Nil] : A6579.6-6585.8 A, [Sii] A6716: 6712.8-6718.1 A and [Sii] A6731: 6727.2-6733.1 A). Fluxes have been 
displayed in a logarithmic scale, in units of 10~^^ erg s~^ cm~^. In all the ma ps, north is up and East is to the left. The HH 110 knots and the spatial scale 
have been labeled in the Ho; panel (30 arcsec ~ 0.07 pc for a distance of 460 pc: lReipurth & 01ber3ll99lh . 




10 ^ 1 10 ^ 1 1 1.5 

Figure 3. Line-ratio maps derived from the integrated line emissions displayed in Fig. [2 The [S Ii]/Ha and [Nii]/Ha maps show the excitation conditions of 
the HH 110 emission, while the [S ii] 6716/6731 map traces the electron density. Contours of the integrated flux of the Ha line have been overlaid to help with 
the knot identification, which have been labeled in the [N Iij/Ho; map. 



Integral Field Spectroscopy ofHH 110 5 




cloud. Contours and knot identification, as in Fig. [3] 



3.1.2 Excitation and density conditions 

First results on the physical conditions of HH 110 were ob- 
tained through long-slit spectroscopy acquired along t he jet outflow 
(iReipurth. Rasa & Heathcotelll 99^jRiera et al"2003a') and at a sev- 
eral positions across the jet beam (Riera et al 2003a). The widths 
of the slits used in these observations were 2 and 1.5 arcsec, re- 
spectively. In both works, the [S ll]/Ha ratios derived for the knots 
intersected by the slit along the jet axis outflow range from ^0.2 
to ~0.7. T hese ratios c orrespond to an intermediate/high degree of 
excitation (iRaga et aL| [l996). A trend was found of an increasing 
[S Il]/Ha ratio (hence a decrease of the excitation) moving from 
knot C to knot M, where the lowest excitation was detected. A de- 
crease of this ratio (hence an increase of the excitation) was found 
beyond knot M. In general, the values derived in both works are 
well consistent, except for knots A and N, where the [S ll]/Ha ratio 
reported by Riera et all (l2003ah is s ignificantlv higher (by a factor of 
two) than the ratio reported bv lReipurth, Raga & Heathcote^ (.199 6). 
This is not surprising since the excitation conditions should change 
significantly within arcsecs across the jet beam, as first revealed by 
t he sp ectroscopy of the jet cross section performed by Riera et al 
(l2003ah . Thus, the sampled regions of knots A and N might not 
coincide in both observations. Furthermore, no clear trends for the 
excitation structur e of the jet cross s ections were found from the 
spectra acquired by'Riera et al' ('2003a'). These authors found a com- 
plex structure at the four positions across the jet beam, where the 
spectra were obtained. 

To look for a more complete picture of the entire spatial 
structure of the excitation and density of the jet, we obtained the 
[Sll](6716-h6731)/Ha and [Nil] (6548-^65 84)/Ha line-ratio maps 
from the IFU data (Fig. O. The line-ratio maps were created by 
dividing spaxel-by-spaxel the corresponding flux-integrated line 
maps (all the spaxels with values below 5 a were blanked before 



to obtain the line-ratio map). At the weakest knots, the estimated 
uncertainties are 10% and 15% in the [S II] (67 16-^6731 )/Ha and 
[Nil] (6548-^65 84)/Ha line-ratio values respectively. 

The [S ll]/Ha and [N ll]/Ha maps clearly reveal the complex 
excitation structure of the emission that was only barely outlined 
from long-slit spectroscopy. Several facts should be remarked: (i) 
for the knots that were sampled through long- slit spectroscopy, the 
line-ratio values derived from the IFS data around the peak posi- 
tions are in good agreement with those derived from the long- slit 
spectra. For knots A and N, the ratios obtained fro m the I FS data 
are in better agreement with those derived bv Riera et al (l2003ah 
than with those derived bv lReipurth. Raga & Heathcot^ l 19961) : (ii) 
IFS line-ratio maps revealed some regions in which the degree of 
gas excitation is low. In particular, the highest [S Il]/Ha ratios (^ 
1.5) are found around knot H and slightly lower ratios, which are 
also compatible with a low degree of excitation, are found around 
knots R and S. This wider range of excitation degree could not be 
detected by previous long-slit observations, which did not sample 
these low excitation regions; (iii) no clear systematic trend in the 
spatial distribution of the excitation was found neither along the jet 
outflow nor across the jet beam. However, the integrated [S ll]/Ha 
and [Nll]/Ha maps indicate that the excitation is higher for the 
northern jet knots A to D (where the outflow has a higher collima- 
tion) than beyond knot E, from the region where the outflow widens 
and the knots begin to lose their alignment relative to the axis di- 
rection traced by knots A-C. For this jet region (i. e. from knots E to 
Q) the line ratio maps show some signatures of an increment in the 
excitation moving outwards, from the knot peak to the eastern knot 
border and to the interknot. The line ratio maps are thus suggestive 
of a decrease in the excitation from east to west across the jet beam. 
This is true for all knots, except around knot I (at the western jet 
border), where the excitation is higher and similar to that derived 
in the northern (A to D) knots. 
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The spatial distribution of the electron density (rie) is traced 
by the [S ll] 6716/6731 map. Figure [3] shows the map created from 
the IFU data (with the assumptions mentioned before for the excita- 
tion maps). At the weakest knots, the estimated uncertainty is 15%. 
The values found range from ^ 0.9 to ^ 1.4, which correspond to 
Tie ranging from 1000 to 50 cm~^ respectively (these values were 
derived using the TEMDEN task of the IRAF/STSDAS package, and 
assuming Te = 10^ K). The electron density cannot be properly 
evaluated around the knot R region, because the derived line ratios 
are higher than the values for which the [Sll] 6716/6731 is den- 
sity sensitive in the low-density limit. In general, the trend found 
is a decrease in the rie, moving along the jet from north to south. 
This result is consistent with the density b ehaviour of the knots 
derived fro m the long - slit spe ctra by Reipurth. Raga & Heathcotel 
(Il996h and lRiera et all (l2003ah . However, the spatial distribution of 
the density obtained through IFS (Fig. [3]) is more complex than that 
outlined from long- slit spectroscopy, and could not be properly de- 
scribed based on this kind of observations. The density map indi- 
cates that, in general, the electron densities are higher around the 
knot peak positions and decrease towards the edges of the knots and 
towards the interknot gas. Across the jet beam, a trend is observed 
in the knots peaking towards the west side of the axis defined by 
knots A-C (e. g. knots E, I, N, P) that are denser than knots peak- 
ing towards the east side of this axis (e. g. knots D, J, O, Q). The 
density map also shows some local departures of this general trend 
that are significant. In particular, it should be noted that the high- 
est Tie is reached around knot E (i. e. ^ 40 arcsec south from knot 
A). In fact, the IFU data revealed an enhancement in the electron 
density, reaching values of ^ 1000 cm~^ and with a maximum 
of 1300 cm ~^ located 3 arcsec southwest of the knot E peak 
intensity. This denser clump, extending over 30 arcsec^, appears 
well delimited at its eastern edge, as shown by a sharp fall of rie, 
which reaches values of 200-300 cm~^ by only moving :^ 2 
arcsec east from the E peak intensity. In addition, although rie gen- 
tly decreases from knots A (-- 600 cm~^) to H (-- 400 cm~^), the 
knot I, southwest of all of these knots, has a rie similar to knot B 
{r^ 500 cm~^), while knot D, southeast of knot C, has similar low 
values (^ 150 cm~^) to those found at the southern, more rarefied 
knots, beyond knot N. 



3.2 Kinematic analysis 

3.2.1 Spatial velocity distribution 

The emission lines at most of the positions mapped in HH 110 
show asymmetric profiles that cannot be properly fitted by a sin- 
gle Gaussian. Hence, we calculated the flux-weighted mean radial 
velocity (or first-order intensity momentum) instead of the line cen- 
troids of a Gaussian fit, and the flux- weighted rms width of the line 
(or second-order intensity momentum), proportional to the FWHM 
of the line, and described the integrated velocity field of the je{3 
from these intensity momenta. HH 110 shows a velocity field be- 
haviour very similar in all the mapped lines (Fig.©. The radial ve- 
locities are blueshifted (Vlsr ranging from '-^ -10 to -80 km s~^ 
in Ha), with an increase towards more blueshifted values mov- 
ing from north to south. Velocities that were derived from Ha are 
blueshifted by ^ 20 km s~ ^relative to those derived from the [N ll] 
and [S II] (the velocities derived from these two lines are similar). 



^ All the velocities in the paper are referred to the local standard of rest 
(LSR) velocity reference frame. A Vlsr = +8.5 km s~^ for the parent 
cloud has been taken from Reipurth & Olberd (19911). 




Figure 5. Maps of the FWHM (in km s~ ^) for the line labeled in the corre- 
sponding panel. Contours and knot identification, as in Fig. [3] Spaxels with 
a FHWM below 115 km s~^have been blanked. (At several positions, like 
in the southern knots, FWHM should be narrower and could not be properly 
determined due to the instrumental line profile). 



The velocity offset found between Ha and [S ll] is not an artifact 
caused by an inaccurate wavelength calibration of the IFS data, 
since it has been also found in long-slit, higher resolution data (e. g. 
iLopez et al .12005). A more detailed inspection of the velocity maps 
reveals a more complex structure of the velocity field: (i) the less 
blueshifted velocities (Vlsr ^ -10 to -15 km s~^ in Ha) are found 
around the northern emissions X, Y and in knot D, all of them lo- 
cated towards the eastern side of the jet axis direction traced by the 
knots A-C; (ii) conversely, more blueshifted velocities (Vlsr ^ - 
50 to -60 km s~^ in Ha) appear around the regions of knot O 
and Q, which also are located at the eastern jet side. Even more 
blueshifted velocities (up to Vlsr ^ -80 km s~^ in Ha) are found 
around knots I and K, located at the western jet border, quite out- 
side of the main jet body. 

The maps of the FWHM for the Ha and [S II] A 6716 A lines 
obtained from the second-order intensity momentum are shown in 
Fig. [5] Maps were corrected from the instrumental profile by sub- 
tracting it in quadrature spaxel-by-spaxel. Spaxels with values be- 
low a 5cr rms level were blanked. The spatial distribution of the 
FWHM shows a complex pattern in both lines. In general, widths 
are broader for Ha than for the [S II] lines. The FWHM of the Ha 
ranges fr om 30 to 90 km s~ ^, being in good agreement with those 
found bv lRiera et 'ai(l2003bh from Fabry-Perot data and with those 
derived bv 'Hartigan et al.l from long- slit, high-resolution 

spectra. The northern knots (A to C) present the largest values (80- 
90 km s~^) of the FWHM, decreasing towards the south of the 
outflow (up to 30-40 km s~^ around knots M and Q). In fact, the 
FWHM around the arc-shaped, southern S emission region should 
be even narrower, but the widths could not be properly evaluated 
due to the instrumental line profile. Like for other jet properties, 
this general trend is not well fulfilled by the whole outflow: a re- 
gion having significant narrow lines, their widths being close to 
detection limit, is also found around knot D, while higher FWHM 
values (90 km s~^, similar to the widths derived around knot A) are 
found around knot I. Concerning the FWHM of the [S ll] lines, their 
values range from 20 to 70 km s~^and are in agreement with the re- 
sults of lHartigan et aLl2009») . The spatial distribution of the FWHM 
values for the [S ll] lines are similar to the FWHM Ha map. Some 
key differences in FWHM between Ha and [S ll] should be men- 



tioned: (i) the widths of the [S ll] hnes drop significantly around 
knot C, where its value diminished by a factor of two relative to 
that found around knot A, while knots C and A have comparable 
FWHMs in Ha; (ii) knots M and Q have very similar FWHM val- 
ues in [S II] and Ha, while the general trend for the rest of the 
knots is the Ha line being wider (by 15-20 km s~^ typically) than 
the [S II] line. 



3.2.2 Peak intensity knot velocities 

In an attempt to compare long- slit and IFS data, we have obtained 
a representative spectrum for each jet knot by extracting the spec- 
trum of the spaxel (aperture of 2 arcsec) that corresponds to the po- 
sition of the knot intensity peak (Fig. (6]). The knots show significant 
differences in their physical conditions, as inferred from their rela- 
tive line intensities. Moreover, the lines present asymmetric profiles 
that in some knots are suggestive of being double-peaked (see, e. g. 
knots H, I) although the profiles cannot be resolved because of the 
spectral resolution of the data. 

Figure [7] displays the radial velocities at the position of the 
knot peak intensity, in both Ha and [S ll] lines, as a function of the 
distance to knot A. The two radial velocity fields are very similar. 
A trend of increasing blueshifted (more negative) velocity values 
with the distance to knot A is outlined. This trend is consistent 
with the velocity behaviour found in previous works (iRiera et all 
l2003bl : iHartigan et all [200^ . The radial velocities of some knots 
present deviations of this general trend. In particular, the isolated 
knot K (found outside the jet body) has a blueshifted radial ve- 
locity more negative than the rest of the knots, and its deviation 
from the general trend is statistically significant. A linear fit of the 
peak radial velocity of the Ha line with distance, including all the 
knots (19), was obtained (correlation coefficient r = 0.61; rms fit 
residual a = 15.1 km s~^). The fit is shown in Fig.[7]by a contin- 
uous straight line, and the la deviation from the fit, by the dashed 
lines. Knot K deviates from the fit by 3a. This deviation is sta- 
tistically significant for a sample of 19 knots. Then, a new linear 
fit was calculated after removing knot K from the sample (i.e. 18 
knots). We obtained values of r = 0.74 and a = 10.5 km s~^. 
The knots that deviate more than 1.5cr from this new fit are knot O 
(1.8a deviation), I (2.2cr deviation), and K shows now a 4.6a de- 
viation. Unfortunately there is no proper motion determination for 
knot K that would allow us to calculate the projection angle of its 
motion. Hence, the current data do not allow us to discern whether 
the kinematic behaviour of knot K, with a radial velocity apprecia- 
bly blueshifted as compared with the rest of the knots, can be only 
attributed to the jet geometry (i.e. to projection effects). Recently, 
iHartigan et aP 12009) imaged the HH 110 region in the H2 2.12 
fim line. They report a weak jet emerging from the source IRSl 
(to the west of HH 110) and extending from north to south. This 
weak jet seems to be pointing at knot K. In conclusion, there are 
some indications (a location isolated from the jet flow, and a kine- 
matics significantly different from the rest of the knots) that knot 
K could be related to IRSl, but its relationship with the HH 110 
outflow cannot be fully discarded based on the results presented in 
this work. 

The full spatial velocity (Vtot) and an upper limit for the an- 
gle between the knot motion and the plane of sky (0=0° when 
the flow moves in the plane of the sky) were calculated from the 
knot pr oper motions, deri ved from multiepoch [S II] narrow-band 
images (iLopez et"aDl2005h and the radial velocities, derived in this 
work from the [S ll] lines at the position of the knot intensity peak 
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Figure 6. Spectra of the HH 110 knots in the spectral windows including 
Ha + [Nil] and [Sii] 6716, 6731 lines. Each spectrum has been extracted 
for the spaxel that corresponds to the peak intensity position (i. e.within an 
aperture of 2 arcsec) and its intensity has been normalized to the Ha peak 
intensity. 



(the typical errors for and Vtot are ^ 6° and ^12 km s~^ re- 
spectively). The results are displayed in Fig. [8] 

The values derived for the inclination angle suggest that the 
flow remains close to the plane of the sky (0 ^ 20° for most of 
the knots). The inclination gently increases from ~ 8° at the 
beginning of the outflow (around knot A) to higher values (^ 30°) 
for distances ^150 arcsec from knot A (around knots R and S, 
coinciding with the region where the outflow has a highly curved 
shape). 

As can be seen from Fig. [8] (lower panel), the spatial ve- 
locity (Vtot) derived at most knots is ^ 100 km s~^. There is a 
trend of decreasing velocity with distance along the outflow, from 
Vtot — 185 km s~^ for knot A to Vtot ^ 75 km s~^ for knot S. 
A few knots present a deviation from this velocity trend that might 
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be significant: knot O, with a higher velocity (Vtot— 190 km s~^) 
than the neighbouring knots (Vtot^ 150 km s~^), and knots E and 
B, having lower Vtot than the neighbouring knots. In particular, the 
lowest spatial velocity is obtained at knot B (Vtot— 25 km s~^). 
This value is far from the spatial velocities of its neighbouring knots 
(Vtot^ 165 km s~^) and even far from the velocities of the slower 
knots. In order to better visualize the velocity behaviour to we are 
referring, a linear fit was calculated after removing knot B from the 
sample (i.e. 14 knots). The fit (correlation coefficient r — 0.53; 
rms fit residual cr=31.8 km s~^) has been drawn in Fig. [8] by a 
continuous straight line, and the la deviation from the fit, by the 
dashed lines. As can be seen from the figure, the spatial velocity 
deviates ^ la for most of the knots. A few knots lie outside the 
region delimited by the dashed lines: knot O (2.2a deviation above 
the fit), knot E (2.3a below the fit) and knot B (4.3a below the fit). 
The region around knot B could correspond to the region where the 
outflow is excavating a channe l throu gh the dense mole cular clump 
(iReipurth. Raga & Heathcotd [19961 : ILopcz et all 12005 ^ and might 
trace the location where the jet collides with the cloud. Hence, sig- 
nificant deviations from the kinematics trend that account for the 
jet/cloud interaction should be expected at this location. This may 
explain the appreciably low spatial velocity found for knot B. Con- 
cerning the other two knots with deviations less significant, knot 
O also has a highly blueshifted radial velocity, (although within a 
la deviation from the radial velocity fit, see Fig. [7]). The low Vtot 
at knot E mostly corresponds to the tangential velocity component 
(see also Fig. [7]) As was mentioned before, the highest values of the 
electron density in the outflow, measured from the integrated [S ll] 
line ratio, are found around knot E location. Thus, similarly to what 
happens around knot B, a stronger interaction between the jet and 
the inhomogeneous ambient material around the knot E location 
might be producing these departures in the kinematics and physical 
conditions relative to those found at the surrounding knots. 

3.23 Line Profile Shapes 

As already mentioned in §3.2.1, the line profiles show some degree 
of asymmetry relative to a single-Gaussian shape that varies from 
knot to knot. This can be seen for the Ha and [S ll] lines of the 
spectra extracted at the position of the peak intensity of each knot 
(Fig. [6]). The analysis of spectral line asymmetries through line bi- 
sectors provides a strong tool to understand the shape and origin 
of line profiles, and is widely used in stellar astr ophysics (see e. g. 
lOueloz et al.l200ll : lDall et al.'2006':'Martmez-Fiorenzano 2008). In 
order to make more clear the asymmetries of the HH 110 line pro- 
files, we calculated the line bisectors for the Ha line at the position 
of the knot peak intensity. The results are shown for each profile 
in Fig. |9] We found knots with some high- velocity emission ex- 
cess, as shown by the tilt of the line bisector towards blueshifted 
wavelengths (knots A, C, D, H, J, L, M, O, P, R and S). We will 
refer to these knots as those having "blue asymmetry". For other 
knots, there is low-velocity emission excess, with the line bisec- 
tor tilted towards redshifted wavelengths (referred as knots having 
"red asymmetry": knots E, I and Q). Finally, a few knots show a 
nearly symmetric profile (knots B and N). Such line-profile be- 
haviour might be explained in part by the combined effect of the 
knot velocity and flow inclination at the knot position. 

Some clues can be o bt ained from the work of 
iHartigan. Ravmond & HartmannI (Il987h . who modeled the 
effect of varying the velocity and the inclination angle of the 
flow in the plane of the sky on the Ha line profiles emerging 
from bow shocks. They found that the asymmetry of the line 
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Figure 7. Radial velocities measured at the position of the peak intensity 
of the HH 110 knots, as a function of distance to the knot A peak intensity 
in Hq;. The black dots correspond to the values derived from Ho; and the 
blue triangles, from the [S ii] lines. The typical value of the velocity error 
is ~ 10 km s~^and the uncertainty in the knot position ~ 1 arcsec, and are 
displayed in the lower-right corner of the figure. The straight lines mark the 
linear fit to the data (solid line) and the 1 a deviation from the fit (dashed 
lines). 
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Figure 8. Total velocities (lower panel) and inclination angle, defined by 
the direction of the knot motion and the plane of the sky (upper panel) as a 
function of distance to the position of the [S ii] intensity peak of knot A. The 
typical values of cf) and Vtot errors are ~ 6° and ~ 12 km s~^, respectively 
and are displayed at the corner of each panel. The uncertainty in the knot 
position (~ 1 arcsec) is smaller than the size of the symbols. Straight lines 
that mark the linear fit to the data (solid line) and the 1 a deviation from the 
fit (dashed lines) are drawn on the velocity panel. 
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profile increases with the velocity of the shock, and with the 
inclination angle. Resolved, double-peaked profiles are only found 
for velocities higher than ~ 150 km s~^and for angles greater 
than :^ 45°. Furthermore, moderate variations on the inclination 
angle result in significant effect in the shape of the line profile. For 
a given angle, the shape of the asymmetry changes from "red" to 
"blue" depending on the velocity: e. g. for low inclination angles 
30 °), a "red asymmetry" is obtained for low velocities 100 
km s~^) and a "blue asymmetry" for higher velocities 200 
kms"^). 

Interestingly, most of the HH 110 knots showing a "blue 
asymmetry" in the line profile correspond to knots having high ve- 
locity (Vtot ^ 160 km s~^, e.g.: A, C, J, O) or inclination angle 
(/) ~ 30° (e.g.: R, S). "Red asymmetry" is observed in knots 
with lower velocities (e.g. E, with Vtot = 80 km s~^and Q, with 
Vtot = 110 km s~^) and lower inclination angles (0 ^ 20°). Note 
in addition that some pairs of knots with close values of Vtot and 
(j) show very similar profiles, as is the case of knots C and J, and of 
knots H and L. 

Hence, there is a general agreement of the Ha line profiles 
observed in HH 110 with the ob tained from the shock models of 
iHartigan. Ravmond & Hartmann I (1987). However, other proper- 
ties of the emission could also being contributing to the line shapes. 
For example, different morphologies of the bow emission, keeping 
all the other parameters constant, can give different line profiles, 
as shown by Schul tz. Burton & Brand (12005 ). They found that the 
bluntness at the apex and the radius of the tail are two important 
factors in determining the line profiles for the case of H 2 lines aris- 
ing from a bow- shock. In fact, Hartiganet ID (l200^ found from 
laboratory jet experiments how the working surface of a deflected 
bow shock present complex emission structures as those observed 
in HH 110. These structures give rise to different asymmetric line 
profiles, depending on the HH 110 knot and hence, on the shock 
velocity and projection angle. 

3.2.4 Kinematics behaviour of physical conditions: channel 
maps 

A more detailed sampling of the gas kinematics was obtained by 
slicing the datacube into a set of velocity channels. Each slice cor- 
responds to a constant wavelength bin of 0.3 A (i. e. a velocity bin 
of- 14kms"^: 13.7 km s"^ for Ha and [Nil], and 13.4 km s"^ 
for [S II]). Hence, we have obtained the channel maps for each of 
the emission lines within the observed wavelength range, namely 
Ha, [N II] A 6584 A and [S ll] A 6716 A , as illustrated in Figs. [TOl 
[TT]and[T2l respectively. Moreover, we also analysed the behaviour 
of the excitation and density as a function of the velocity. 

The emission detected above a 3cr SNR level spreads over 
different velocity ranges, depending on the line. Emission in Ha 
was detected over a wider velocity range (:^ 200 km s~^, from 
-135 to +15 km s"^) than in [Sll] (^ 180 km s"\ from -115 
to +10 km s"^) and in [Nil] (^ 150 km s"\ from -100 to 
-1-50 km s~^). Note, however that the morphology of the emission 
does not significantly differ among the three lines when channels 
of similar velocities are compared. 

Ha emission covering the whole velocity range was only de- 
tected at the northern, more collimated flow region, from knots A 
to H. For velocities ranging from -80 to -10 km s~^, there is 
emission detection from all the knots, the stronger emission com- 
ing from the :^ -50 to -23 km subchannels. Emission from veloc- 
ities more blueshifted than ^ -90 km s^^or more redshifted than 
^ -1-20 km s^^were not detected (up to a 3 a snr level) beyond knot 




Relative Wavelength (A) 

Figure 9. Ho; line profiles of the knots labeled in each panel, obtained from 
the spectra extracted at the position of the knot peak intensity (see Fig.|6j. 
The line bisector has been drawn for each profile. In each panel, the line 
intensity has been normalized to the value at the peak. On the x-axis, the 
wavelengths are relative to the rest Ho; wavelenght. 



Q, from the loci where the flow trajectory appreciably changes its 
path. The kinematic behavior in the [N II] emission is similar to that 
found in Ha, although the velocity range with emission detection is 
different. In this sense, the [N II] emission arising from most of the 
knots covers a range from ~ -60 to -I-IO km s"^, the stronger emis- 
sion coming from the -45 to -18 km subchannels. [Nil] emis- 
sion was not detected from the southern knots, beyond knot Q, with 
the exception of some hint (^ 2cr snr level) at -18 km s"^ around 
knot R. Note in addition the lack of [N ll] emission from knots A 
to C for velocities more blueshifted than c^i -90 km s~^. Finally, 
we detected [S ll] emission from the whole flow covering a range 
from -90 to -1-I8 km s~^, the stronger emission coming from the 
:^ -35 to -8 km subchannels. As in the case of Ha and [N ll] , there 
are neither emission more blueshifted than -90 km s^^nor more 
redshifted than ~ -1-20 km s"^ detected beyond knot Q. 

Line-ratio channel maps were created from the IFU data (as 
mentioned in §3.1.2 for the integrated line-ratio maps) to anal- 
yse the excitation and density as a function of the velocity. The 
estimated uncertainties are better than 25% in the [Nll]6584/Ha 
and 10-15% in the [S ll] (67 16-^6731 )/Ha and [Sll] 6716/6731 
maps. Figures [13] and [14] display the [S ll](6716-h6731)/Ha and 
[Nll]6584/Ha line-ratio channel maps, tracing the gas excitation 
as a function of the radial velocity. 

The line-ratio maps for each velocity channel follow the spa- 
tial behaviour of the excitation that was already found in the inte- 
grated line-ratio maps (Fig. O. The gas excitation is higher (i. e. 
lower [Sll]/Ha line-ratio values) at the northern knots (A-D), 
and becomes lower (i. e. the [S ll]/Ha line-ratio values increase) 
beyond knot E, along the region where the flow cross section 
widens. Furthermore, the trend found for this region (a decreas- 
ing excitation moving across the jet beam, from the eastern jet 
side to the west) is now more clearly seen on each of the chan- 
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Figure 10. Ha channel maps centred on the Vlsr radial velocities 
(km s~^) labeled in each panel. The spatial scale is given in one of the 
panels. Fluxes have been displayed in a logarithmic scale, in units of 
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nel map. An opposite trend of increasing excitation from east 
to west is outlined at the curved, terminal S region from the - 
45 km s~^ <Vlsr<+5 km subchannels, where the signal is large 
enough to evaluate the ratio with confidence. 

In addition, the line-ratio channel maps indicate that the gas 
excitation varies depending on the gas kinematics, namely higher 
excitation (i. e. lower [S Il]/Ha line ratios) for the higher velocities 
(in absolute value, both blueshifted and redshifted) and lower ex- 
citation (higher line ratios) for lower (in absolute value) velocities. 
The differences are found for most of the knots. For example, for 
knot C, at the position of its peak intensity, we measured a [S ll]/Ha 
ratio ^ 0.2 for VisR ^ -90 km s"^ and VisR ^ +30 km s"\ while 
a ratio of ^ 0.5 is found for -50 km s"^ < Vlsr < +20 km s~\ 
The same behaviour appears for low-excited knots: at the peak 
intensity position of knot M, the [Sll]/Ha ratio is ^ 1.1 for - 
30 km s"b < VisR < +20 km s"\ being ^ 0.5 for VisR ^ - 
90 km s"b and VisR ^ +30 km s"\ 
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Figure 11. Same as Fig.[Tol but for the [N ii] 6584 A emission line. 
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Figure 12. Same as Fig.[To] but for the [S ii] 6716 A emission line. 
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Figure 13. The [S ii]/Hq; line-ration maps, tracing the gas excitation, de- 
rived from the channel maps at the Vlsr labeled in the corresponding 
panel. 



The [Six] 6716/6730 line-ratio channel maps (Fig. (TS]) trace 
the electron density of the emitting gas as a function of the 
radial velocity. In general, the spatial distribution of the den- 
sity follows the same trends that were found from the integrated 
[S II] 6716/6730 line-ratio (Fig.O. In this sense, density decreases 
from north to south along the jet, but with some exception, e. g. 
around knot E, where density is the largest. Variations on the 
electron density also depend on the gas kinematics. For a given 
knot, the largest electron density corresponds to the emission mov- 
ing at lower (in absolute value) velocities (typically within the 
range -25 km s~^<Vlsr<+5 km s~^), whereas the density falls 
down to 40-50% for the emission moving at high velocity 
(Vlsr^ +20 km s~^and Flsr^ -60 km s~^). Hence, the emis- 
sion from slower gas appears denser and less excited than the faster 
(red or blueshifted) emission, which is more excited and rarefied. 
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Figure 14. Same as Fig. [13] but for the [N ii]/Hq; line-ration maps. 



4 SUMMARY AND CONCLUSIONS 

To conclude, the IFS data confirm the peculiar nature of HH 110 
as compared with most of the stellar jets already observed. The 
most relevant results presented in this paper might be summarized 
as follows: 

• These IFS data have allowed us to generate the first [Nil] 
narrow-band image of HH 110 currently obtained, which shows 
a similar jet cross section to the [S II] emission. In addition, these 
data allow us to confirm that the faint filamentary emission labeled 
X and Y, detected to the east of knots A-C, which probably cor- 
responds to a "fossil" outflow channel, mainly should arise from 
"true" Ha line emission. The contributions from nearby continuum 
or from the [N II] lines, if any, should be very weak, as no emission 
was detected in [Nil], and the continuum contribution was previ- 
ously removed. 

• The line-ratio maps derived from the integrated line intensi- 
ties revealed complex spatial structures for the gas excitation and 
electron density. Although the line ratios that trace gas excitation 
are indicative of an intermediate/high-excitation degree for most of 
the emission mapped, several regions of low-excitation emission 
are detected through the flow (e.g. [Sll]/Ha ^1.5 around knot 
H). In general, the [S ll] line ratios correspond to low electron den- 
sity values (rie ^ 1000 cm~^) and show smooth spatial variations. 
However, strong local enhancements of rie relative to its surround- 
ings were detected at several locations (e. g. around knot E, where 
the Tie value is about three times that of the surroundings). 

• Radial velocities appear blueshifted, relative to the Vlsr of 
the cloud, the modulus increasing with the distance to knot A, 
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Figure 15. The [S ii]67 16/6731 line-ratio maps, tracing the electron density, 
derived from the channel maps at the Vlsr labeled in the corresponding 
panel. 



• The channel maps showed that the emission spreads over dif- 
ferent velocity ranges, depending on the line considered. In spite 
of that, the maps also showed that the morphology of the emis- 
sion in each of the velocity channels is similar for all the lines, 
the strongest emission corresponds to the velocity interval ^- 
20<yLSR<+5 km s~^, where significant emission from all the 
knots was detected. In addition, the previously unknown kinematics 
of the excitation and electron density along the full spatial extent of 
the outflow were obtained from the channel maps. The degree of 
excitation and the electron density appear to be anticorrelated with 
the velocity modulus, in the way that the emission from the gas 
moving faster (either, at red or blueshifted velocities) has higher 
excitation and lower density as compared to the emission from the 
gas moving at slower velocities, which appear denser and less ex- 
cited. 

• From the IFS spectra extracted at appropriate spaxels, we were 
able to compare the physical conditions obtained from these data 
with those derived from previous long- slit and Fabry-Perot obser- 
vations. The general trends followed by the properties of the emis- 
sion and even the derived values of the observables are in good 
agreement at all the commonly sampled positions. However, the 
complex structure of the HH 110 kinematics and physical con- 
ditions, already outlined in previous works, has been now better 
characterized from this more complete IFS sampling. These data 
revealed significant variations of the kinematics and physical con- 
ditions over short distances that were not properly sampled by the 
long- slit data. 

• The emission lines have asymmetric profiles. The analysis of 
the Ha line profiles performed through the line bisectors method 
shows that the departure of the profile from a single-Gaussian shape 
varies from knot to knot. We found that the shape of the profiles 
seems to be related to the full velocity and the viewing angle of the 
knot. 



although with some appreciable local departures from this trend 
(in particular, knot K). This systematic increase in radial velocity 
should not be interpreted as a sign of jet acceleration, since no simi- 
lar trend was derived for proper motions. In fact, the trend found for 
tangential velocities goes in the opposite way (decreasing velocity 
as a function of distance from knot A, with some appreciable local 
departures, e. g. in knot B). Thus, this apparent acceleration should 
be better attributed to a geometric effect: the interaction between 
the jet and the inhomogeneous ambient surroundings will lead to 
the deflection of the jet and therefore, to changes in the projection 
of the full spatial velocity along the line of sight (radial velocity) 
and along the plane of the sky (tangential velocity). 

• In general, the line profiles are broader for Ha than for [S II], 
and the FWHM in both lines decreases from the northern to the 
southern knots, although this trend is not fully systematic along the 
outflow. 

• We recalculated the full velocity and the motion direction rel- 
ative to the plain of sky of the HH 110 knots using the [Sll] ra- 
dial velocities, derived from the spectra extracted at the position 
of the knots peak intensity, together with knot proper motions, ob- 
tained from ground-based, multi-epoch [S ll] narrow-band images. 
We confirmed the general trend of a decrease in the velocity and 
an increment in the projection angle with distance from knot A, as 
found in previous works that used a more limited knot sample. Sig- 
nificant departures from the trend are found in some knots. These 
knots could be tracing the loci where a stronger interaction between 
the jet and the inhomogeneous medium is taking place. 



The scenario first outlined by iReipurth. Rasa & Heathcotel 

(Il996h (i. e. deflection of t he HH 270 j et by a dense molecular 
cloud) and later modeled bv lRaga etaLl(l2002h sives a reasonable 
explanation on the origin of the HH 1 10 jet and offers a good qual- 
itative fitting for a set of properties observed in the jet behaviour 
(e. g. the strong lack of symmetry, relative to the outflow axis, of 
the proper motions), but fail in reproducing in more detail the kine- 
matics and some other prope rties through the outflow. In particular, 
as was already discussed b 
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Lopez et aD |2005), the jet/cloud col- 
2002') did not predict the deceleration 
of the full spatial velocity along the outflow that is derived from 
imaging, long- slit and IFS data. 

A variability in the ejection outflow velocity (i. e. pulsed jet 
models) and/or an anomalously strong interaction between the out- 
flow and the inhomogeneous environment were then suggested as 
plausible alternatives giving rise to such kinematic behaviour. In- 
terestingly, new modeling that incorporate s the se mechanisms has 
been recently published. Yirak et al.l (l2008h and lYirak et aP (l2009h 
performed more complex, two-dimensional hydrodynamical simu- 
lations of jets that propagate through a small-scale inhomogeneous 
environment (e. g. clumps with a size smaller than the jet beam). 
Their simulations predict that the collisions between slow clumps 
overtaken by faster ones should produce shock structures slightly 
smaller than the jet beam and displaced from the jet axis, like what 
is observed at some locations along HH 110. These simulations 
also predict a complex evolution of these shock structures. In par- 
ticular, properties such as excitation and velocity profiles across the 
jet beam, are far more complex than expected from an unperturbed 
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clump, or from the internal working surfaces arising from current 
models of a pulsed jet. 

Further insights on the beha viour of d eflected supersonic jets 
were recently obtained by Hart igan et al ] (|20Q9) from laboratory 
experiments. The experiments show how the morphology of the 
contact discontinuity between the jet and the obstacle develops a 
complex structure of cavities that gives rise to the entrainment of 
clumps of material from the obstacle into the flow. The experiments 
also succeeded in reproducing the morphology of some filaments 
such as those observed in high- spatial resolution (HST) images 
of HH 110, although the observed dynamics in these filamentary 
structures were not properly reproduced by the experiments. Ac- 
cording to the results of these laboratory experiments, as compared 
to their high-resolution, long- slit data, the authors proposed that 
the best model for HH 110 still is that of a pulsed (from a varying 
driving source) jet interacting with a dense molecular clump. The 
results derived from our IFS observations that can be compared 
with those found from these authors (e. g. the velocity behaviour 
through the flow and the line profiles at selected positions) appear 
well consistent, thus giving further support to models that involve 
deflection of a pulsed jet propagating through an inhomogeneous 
ambient medium. 

The IFS data here presented give a full spatial coverage in the 
Ha, [N II] and [S ll] emission lines of the singular outflow HH 110. 
They have allowed us to explore for the first time the whole spa- 
tial distribution of the physical conditions and its relationship with 
the kinematics of the jet emission. We would like to point out that 
there are very few IFS observations of stellar jets with a wide spatial 
coverage of the o utflow in several emission lines (as e. g. HH 34 by 
Beck et al .120071) . Such observations are highly necessary to under- 
stand the behaviour of the physical conditions as a function of the 
kinematics, as well as to explore whether this behaviour varies (as 
in HH 34) or not (as in HH 1) with th e distance from the exciting 
source (see iGarcia-Lopez etal Unfortunately, the paucity 

of available data prevents to establish a general picture on these 
topics. 

Because of the rather chaotic behaviour of HH 110 outlined in 
previous works, as derived from more limited narrow-band imag- 
ing and long-slit spectroscopic data, the IFS observations discussed 
here are particularly useful for characterizing the properties of the 
whole outflow. Hence, a more realistic picture has arisen, suitable 
for designing new state-of-the-art simulations to match the HH 110 
scenario. Note that IFS data give much more detailed information 
on the spatial distribution of excitation as a function of velocity 
than provided by current model simulations or laboratory jet ex- 
periments. These data could provide valuable clues to constrain the 
space parameters in future theoretical works, which are necessary 
to understand the origin, structure and dynamics of HH 110. 
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